QCD FITS TO DIFFRACTION * 
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The most recent HI and ZEUS data for diffractive structure functions are analyzed under three 
different theoretical approaches. This includes the Pomeron Structure Function (PSF) framework, 
Bartels-Ellis-Kowalski-Wiisthoff (BEKW) color dipole approach and the Golec-Biernat-Wiisthoff 
(GBW) saturation model. The models are shown to successfully fit the set of combined data. 
Conceptual differences between the models are discussed and, as an example, the prediction for 
the longitudinal diffractive structure function is presented. 
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1. Introduction 

The process of inelastic electron-proton 
diffraction, ep — > eXY, have been widely 
studied over last years at the HERA collider. 
Quite recently both HI and ZEUS collab- 
orations published particularly precise data 
[1-4]. In this short letter we present the re- 
sults of the analysis of this data in the frame- 
works of three different theoretical models. 
These are: the Pomeron Structure Func- 
tion (PSF) model formulated in the frame- 
work of Regge phenomenology, the Bartels- 
Ellis-Kowalski-Wusthoff (BEKW) two gluon 
exchange dipole model and the satura- 
tion model of Golec-Biernat and Wusthoff 
(GBW). By fitting each of these models to 
the combined set of all measurements we test 
also the compatibility of the data obtained 
with various experimental methods. 

When we plot the experimental data for 
the diffractive structure function [1-4] we ob- 
serve that it decreases in the range of low xp 
values, reaches a minimum at certain point, 
and starts to grow up as xp goes to 1. Thus, 
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F®^ (Q 2 ,P,xp) is expressed as a sum of two 
components 
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In our analysis we concentrate on the low-xp 
part of the diffractive structure function and 
test various theoretical frameworks in which 
this quantity may be calculated. The large- 
xp part is modeled as a Reggeon exchange 
with the structure function of the pion [5] 
and only its normalization for each data set 
is determined from the fit. 

2. Data sets and fitting method 

We use the data for the reduced cross sec- 
tion a^ 3 \Q 2 , /3, xp) or the diffractive struc- 
ture function F^ 3 \q 2 , (3, xp) published by 
HI [1, 2] and ZEUS [3, 4] collaborations, 
respectively. These data sets are, however, 
obtained with different methods and there- 
fore for different upper limits on the proton 
dissociation system mass My. Hence, they 
need to be normalized to the same experi- 
mental situation, My < 1.6 GeV, in accor- 
dance with [1] . Taking all the above into ac- 
count we have the following data sets which 
are used in the fits 
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Table 1. Results, in terms of x 2 > of the fits to the three data sets for the PSF, BEKW 
and GBW models. 



X 2 /(nb data points) 


H1RAP 




ZEUSMX 


All data sets 


PSF 


250.3/240 = 


1.04 


100.9/102 = 0.98 


377.6/444 = 0.85 


BEKW 


286.8/247 = 


1.16 


193.5/142 = 1.36 


493.1/493 = 1.00 


GBW 


272.9/247 = 


1.10 


268.0/142 = 1.88 


564.5/493 = 1.15 



called H1RAP, default data set, not cor- 
rected further, 

• Frf^ measured by the ZEUS using the 
M x method [3] with M Y < 2.3 GeV called 
ZEUSMX, multiplied by the factor 0.85, 

• ay^ 3 ' measured by the HI with the proton 
detected in roman pot detectors [2] called 
H1TAG, multiplied by the factor 1.23, 

• measured by the ZEUS with the pro- 
ton detected in roman pot detectors [4] 
called ZEUSTAG, multiplied by the fac- 
tor 1.23. 

We perform three fits for each model. 
The fit to HI RAP data alone, the fit to 
ZUESMX data alone and the combined fit 
to all four data sets listed above. In the first 
two cases only statistical and uncorrelated 
systematic uncertainties added in quadrature 
were used in the calculation of x 2 whereas 
the total errors were taken in the third case. 
The cut Q 2 > 4.5 GeV 2 was imposed for all 
fits. 

3. Pomeron Structure Function 
(PSF) model 

In the framework of Regge phenomenol- 
ogy the low-xp component of the diffrac- 
tive structure function is attributed to the 
Pomeron exchange. The ip is factorized 
from (3 and Q 2 dependence and we have 

^(3),low-x P=/p/p(a;]p) ^(P) (g2)/3))(2) 

where the first factor is called the Pomeron 
flux. It has the standard form [6] with the 
Pomeron intercept ctp(O), which governs the 
xp dependence and is let as a free parameter 



in our fit. The second factor F 2 {Q 2 , 0) is 
interpreted as the Pomeron structure func- 
tion defined in analogy to the proton struc- 
ture function. In the same manner diffrac- 
tive parton distribution functions (DPDFs) 
arc defined for the Pomeron. The evolution 
of these distributions with Q 2 is governed by 
the DGLAP [7] equations. 

In our analysis wc take the following 
form of singlet zS and gluon zG distributions 
at initial scale Qq 

zS(z,Ql)=A s z Bs (l-z) Cs 

■{l + D s z + E S s/z~)-e^, 
zG(z,Q 2 ) = A G (l-zf°-e^. (3) 

They are similar to that used in fit A of HI [1] 
however we increased the number of param- 
eters adding Ds and Es- This allowed us 
to make the fit independent on the choice of 
Qq, which is fixed at Q 2 , = 3.0 GeV 2 . Hence, 
we have eight parameters introduced by the 
function defined in Eq. (2). The charm con- 
tribution is not fitted but calculated accord- 
ing to [8]. Since in this approach the higher 
twist effects are not present we need to re- 
strict the kinematic range by introducing two 
additional cuts M x = 2.0 GeV and < 0.8. 

The fit results in terms of x 2 are pre- 
sented in the Table 1. Very good fit quality 
is obtained in all three cases. The gluon den- 
sity turns out to dominate in the Pomeron for 
all fits however its shape varies between fits 
with different data sets. We obtain ap(0) = 
1.118±0.005 in the combined fit which agrees 
with the result found by HI (fit A) [1]. The 
full set of parameters as well as the singlet 
and gluon distributions can be found in [6]. 
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4. Bartels-Ellis-Kowalski- 
Wusthoff (BEKW) model 

In this model proposed in [9] the photon- 
proton interaction is realized in two stages. 
In the first stage the virtual photon splits 
into qq or qqg Fock states known as color 
dipoles. Then the dipole interacts with the 
proton via two gluon exchange. Thus, the 
simplest model of Pomeron as the gluon 
pair is assumed. The low-xp component of 
the diffractive structure function consists of 
three terms 



pU W , Iow-xp = pgg + pgg + pqqg ^ 



where T (L) refers to transverse (longitu- 
dinal) polarization of the incoming photon. 
From the properties of the wave function the 
(3 dependence of each term was deduced. The 
qq~T term dominates for intermediate (3 val- 
ues whereas qq~L is the most important at 
large (3. Moreover, the longitudinal compo- 
nent behaves like a higher twist part so nei- 
ther (3 nor Mx cuts are needed. For the 
low (3 values the qqgr term is the largest. 
The xp dependence is impossible to guess 
from perturbative QCD therefore it is as- 
sumed in the form {xo/xp) ri2ii '>^ ) where 
n 2 (4)(Q 2 ) = n 2{i) + n\ {i) \n[Q 2 /Q 2 + 1] with 
n 2(4)' n 2(4) to be determined from the fit. 
Altogether the function given in Eq. (4) has 
eight parameters. 

In the second row of Table 1 we present 
the x 2 values resulting from the fit of the 
BEKW model to the three data sets. Good 
fit quality is observed in the case of HI RAP 
and the combined fit. The fit to the 
ZEUSMX data set gives slightly worse x 2 - 
For all details concerning the parameter val- 
ues and plots illustrating (3 dependence of 

p D(3), low-zp . n BEKW modd the reader 

is referred to [6]. 



5. Golec-Biernat-Wiisthoff 
(GBW) saturation model 

The saturation model proposed in [10, 11] 
was formulated in the color dipole picture. 
Hence, the general structure of F®^' low_xlp 
is identical with the BEKW model and there- 
fore given by Eq. (4) . However in this model 
the cross section for the interaction of the 
color dipole with the proton is assumed in 
the form 

a(x P ,r) = <T {l - cxp(-r 2 Q 2 (x P )/4)} , 

(5) 

where Q 2 s {xp) — (xo/xp) is called the sat- 
uration scale. So, the cross section satu- 
rates either for r — > oo and xp fixed or 
xp — > and r fixed. We used the version of 
the GBW model with three massless quarks. 
The model has only four parameters to be 
determined from the fit (see [6]). No (3 or 
Mx cuts are imposed for the same reason as 
discussed in the previous section. 

We present the results of x 2 from the 
GBW fits in the third row of Table 1. Good 
quality is obtained in the case of H1RAP and 
the combined fit. Significantly worse result is 
found when the GBW model is fitted to the 
ZEUSMX data set. The results for parame- 
ters, especially xq, vary significantly depend- 
ing on the data set used in the fit. Again, we 
refer the read to [6] for further details. 

6. Conclusions 

We have presented the results of the fits of 
three different theoretical models to the most 
recent data on the diffractive structure func- 
tions from HI and ZEUS collaborations. We 
have shown that the combined fit to the four 
data sets with total errors is successful, in 
terms of x 2 , for all models. Especially good 
description of the diffractive data is found in 
the framework of the PSF and BEKW mod- 
els. In addition, the PSF approach works 
good also when fitted to the ZEUSMX data 
set, which is not the case for BEKW and 
GBW models. 
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Let us stress that the three frameworks 
presented here are based on significantly dif- 
ferent theoretical concepts. This is reflected 
e.g. in the twist structure of the models. As 
we see in Fig. 1, the predictions for the longi- 
tudinal part of the diffractive structure func- 
tion Ff are very different between PSF and 
BEKW/GBW approaches. Since in the PSF 
approach only the twist-two component is 
present we expect it to give correct Ft in the 
low (3 region. On the other hand, because Ft 
in the BEKW/GBW models has only twist- 
four part this models are reliable for large f3. 

Finally, let us mention that in [6] the 
reader may find also the results for the 
Bialas-Peschanski (BP) model [12] based on 
the BFKL Pomeron approach. The com- 
bined fit of this model in the kinematic 
range as in the case of the PSF model and 
additional cut Q 2 < 120 GeV 2 leads to 
X 2 /(nb data points) = 1.26. 

For more details concerning the analysis 
described here we refer to the original pa- 
per [6]. 
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Fig. 1. Predictions for the longitudinal part of the 
diffractive structure function as a function of /3. 
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